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1. Introduction 

The spectroscopy of excited hadrons still poses a significant problem to the lattice QCD com- 
munity. Since a Euclidean-time hadron correlator is composed of a sum of contributions from the 
different mass eigenstates, each of which being exponentially suppressed by a factor of the mass, 
disentangling the higher mass states from those of the dominant ground state (let alone the other 
excited masses) can be a daunting task. 

In the following, we present a method (see also [|I|]) for handling just this "disentanglement" 
of masses and we show our latest results for the meson sector. In another contribution to this 
conference [||], we present the results for baryons. 



2. Method 

Our method is based upon the variational approach developed by Michael [||] and later refined 
by Liischer and Wolff [Qj. We use a number of source and sink interpolators to create a matrix of 
correlators which we then diagonalize. 

The basic structure of our meson interpolators is 

O — u nw T* d nw , (2.1) 

where we use T = 75 and 7475, 7 and 747, 757, 77,, and 1 for pseudoscalar (PS), vector (V), 
pseudovector (PV), tensor (T), and scalar (SC) mesons, respectively. The subscripts n and w refer 
to narrow and wide quark sources. These are created using gauge-covariant Jacobi smearing to 
approximate Gaussian distributions of widths 0.27 and 0.41 fm. By opening this new degree of 
freedom to the mesons, we hope to improve overlap with states exhibiting radial excitations (a 
difference in sign between the narrow and wide source contributions will create a radial node in 
the corresponding quark wavef unction). We use degenerate u and d quark masses. For strange 
mesons, we replace one of the light quarks by the strange quark (the physical strange quark mass 
is determined via the PS kaon). 

For our non-strange PS and V mesons we then have up to 6 different interpolators (nw = wn for 
degenerate quark masses) for our basis, giving up to a 6 x 6 correlator matrix (this becomes 8x8 
for the strange mesons). For the PV, T, and SC mesons, we have a basis of up to 3 interpolators (4 
for the strange mesons). 

After creating our correlator matrix, C(t), we must then solve the generalized eigenvalue prob- 
lem: 

C(t)v^ = A«(OC(f )v W , (2.2) 

where the X^ k '(t) are the eigenvalues corresponding to the eigenvectors yW. We note in passing 
that we arrive at almost identical eigenvalues, up to a constant factor, with a normal eigenvalue 
problem, i.e., by omitting the matrix C(to) on the r.h.s. (the eigenvectors, however, are a different 
matter). The eigenvalues provide the time dependence, and hence the mass, of the separate mass 
eigenstates: 

oc e - tM '[\+tf{e- ,AMk )} . (2.3) 

The eigenvectors provide the corresponding wavefunctions (within the highly limited basis). The 
0'(e~ tAMk ) terms in Eq. (2.3) can be reduced by changing, even by further limiting, the starting 
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Figure 1: Effective mass plots for pseudoscalar and vector mesons from our coarse (16 3 x 32, with am a — 
0.05,0.1,0.2 from top to bottom) and fine (20 3 x 32, with am q = 0.04,0.08,0.16) lattices. Both ground and 
excited states are shown, along with the M ± a u results (horizontal lines) from fits of the eigenvalues in the 
corresponding time intervals. 



basis. We observe such effects in our results and, in the end, see the best effective mass plateaus 
in the excited-state eigenvalues by choosing 3x3 correlator matrices throughout (nFn, nFw, wTw; 
the only exceptions are for the PS and V mesons on our coarse lattice, where we use 4 basis states: 
nTn, ny^Tn, ny^Tw, wy^Tw). 

We have results from two sets of quenched configurations, created using the Liischer-Weisz 
gauge action. The lattice spacings are determined via the Sommer parameter [§]: a = 0. 148 fm for 
the 16 3 x 32 and a = 0. 119 fm for the 20 3 x 32 (both spatial volumes are rs 2.4 fm). The quark 
propagators are created using the chirally improved (CI) action [0]. 



3. Results 

In Fig. [T[ we show some of our effective mass plots for the PS and V mesons. Results are 
shown for both the first and second eigenvalues, corresponding to the ground and first-excited 
states. From the apparent plateaus, we choose appropriate time intervals over which we fit the 
eigenvalues with an exponential and extract the mass. 
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Figure 2: Excited-state pseudoscalar masses vs M\ for both lattice spacings. The quark masses are degener- 
ate. All quantities are in units of the Sommer parameter r$. The diamond represents the experimental point 
and the vertical lines mark the values of r\M\ corresponding to the physical strange quark mass. 



Repeating this excercise for many quark masses, we arrive at plots like that of Fig. ||, where 
the excited-state PS mass is displayed as a function of the ground-state PS mass squared (oc m q ). 
Although the statistical errors are substantial, the masses appear to be extrapolating to a value 
consistent with the physical one. Also, the results from the two lattice spacings appear to be 
consistent with each other, suggesting small discretization effects for this state. 

In Fig. |H similar plots are shown for the V, PV, T, and SC states. For the V mesons, the results 
are roughly consistent with the experimental values, although it is apparent that discretization ef- 
fects are larger for these excited states. The other states are more problematic. The excited PV is 
obviously higher than one might naively expect and results for the other excited states are virtually 
non-existent (few reliable effective mass plateaus could be found). We plan to improve the picture 
for such states by including more "realistic" interpolators (e.g., p-wave sources). For the SC meson 
at small quark masses, ghost states dominate the correlation function (which then becomes negative 
[^]), and we therefore plot only the results at large quark masses. 

Figure ^ shows the same plots for the strange PS and V mesons. Again, the excited PS results 
are consistent with the experimental value and display small discretization effects. The excited V 
mesons, however, appear to remain high and also suffer more from discretization effects. 

When looking at the eigenvectors in the 3x3 smearing basis alone, invariably we find the 
same sign for all components in the ground state but sign differences for the excited state. In other 
words, within this limited basis, the excited state has a radial node in its wavefunction. Including 
different spin structures (r and /4r) in the PS or the V basis, the interpretation becomes more 
ambiguous. 

Of course, all of these results are from quenched lattices with relatively small volumes (« 2.4 
fm). So we have some systematic effects which are rather difficult to quantify without repeating 
our analysis on larger, dynamical configurations. It is quite possible that the small volume may be 
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Figure 3: Ground- and excited-state meson masses vs M\ for both lattice spacings. Vector, pseudovector, 
tensor, and scalar mesons appear. The quark masses are degenerate. All quantities are in units of the Sommer 
parameter r§. The diamonds represent the experimental points and the vertical lines mark the values of r\M\ 
corresponding to the physical strange quark mass. 



the reason the excited V and PV mesons remain higher than the experimental values. 
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Figure 4: Ground- and excited-state meson masses vs for both lattice spacings. Both pseudoscalar and 
vector mesons are shown. One quark mass is fixed to the physical strange quark mass. All quantities are 
in units of the Sommer parameter ro- The diamonds represent the experimental points and the vertical lines 
mark the values of r^M\ corresponding to the physical strange quark mass. 
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